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ABSTRACT 

This  thesis  examines  the  influence  of  target  location 
error  upon  small  arms  weapons  system  evaluation.   The 
adequacy  of  the  diffuse  target  approximation  is  examined 
by  comparison  with  tabulated  results  for  a  salvo  of  N-rounds 
of  small  arms  fire. 
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SYMBOL  DEFINITION 


a  Lethal  radius  of  fragmentation 

projectile 

f (XD,Y  *x,  ,y.  )  Probability  density  function 

of  ballistic  error 

f(X.,Y.;x  ,y  )  Probability  density  function 

of  aiming  error 

f(XT,Y..;x  ,y  )  Probability  density  function 

of  target  location  error 

£(XI,YI;XB,YB)  Lethality  function 

N  Number  of  projectiles 


PR(N)  Salvo  kill  probabilit 
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I.   INTRODUCTION 

A.  BACKGROUND 

The  trend  in  small  arms  development  has  been  away  from 
cumbersome,  large  caliber,  slow  firing  weapons  to  light- 
weight, small  caliber,  rapid  firing  weapons.   Sir  Basil 
Liddell  Hart  [18]  stresses  the  roles  played  by  deception 
and  mobility  in  all  warfare,  both  ancient  and  modern.   The 
decisive  roles  played  by  deception  and  mobility  have  made 
target  detection  (location)  an  extremely  important  variable 
to  consider  when  evaluating  and  hence  selecting  a  small 
arms  weapon  system. 

In  this  thesis  the  influence  of  target  location  error 
upon  several  common  measures  of  effectiveness  (MOE's)  used 
in  small  arms  weapon  system  evaluation  is  examined.   The 
adequacy  of  the  diffuse  target  approximation  [7]  is  also 
examined  by  comparison  with  tabulated  results  for  salvo- 
firing  of  N-rounds  of  small  arms  ammunition  against  a 
square  target. 

B.  MEASURES  OP  EFFECTIVENESS  (MOE) 

In  view  of  the  large  number  of  small  arms  weapon  systems 
available  to  the  armed  services  today,  an  acceptable  and 
practical  set  of  criteria  for  comparing  the  overall  combat 
effectiveness  of  these  various  systems  is  of  great  importance 
The  following  measures  were  approached  in  this  thesis:  (1) 
hit  probability,  (2)  lethality,  (3)  kill  probability,  and 


(k)    rate  of  fire.   These  measures  were  selected  because  of 
their  wide  usage  and  ease  of  representation  by  analytical 
models.   References  3  and  13  are  recommended  to  those 
interested  in  additional  measures  sometimes  used.   These 
references  also  contain  a  short  discussion  of  the  problems 
encountered  when  trying  to  decide  which  measures  are 
important . 

C .  PURPOSE 

The  purpose  of  this  thesis  is  to  develop  a  means  of 
effectively  using  the  measures  listed  earlier  when  comparing 
various  small  arms  weapon  systems. 

D.  APPROACH 

The  approach  taken  in  this  thesis  is  to  incorporate 
target  location  error  into  the  various  analytical  expressions 
used  with  the  measures  listed  earlier.   The  effects  of  loca- 
tion error  upon  these  measures  is  th< '  ted  to  obtain 

meaningful  answers  to  the  following  questions: 

1.  When  is  it  desirable  to  sacrifice  accuracy  for  an 
increase  in  the  rate  of  fire? 

2.  When  is  it  desirable  to  use  ball  ammunition, 
shotgun  ammunition,  or  grenades? 

3.  Which  system  best  incapacitates  the  target? 

4.  Which  systems  provide  the  best  suppressive  fires? 
Target  location  error  was  selected  to  permit  examination 

of  the  survivability  of  the  individual  foot  soldier  under 
both  static  and  dynamic  conditions.   What  type  weapon  system 


is  most  desirable  when  very  little  is  known  about  the  true 
target's  location?  What  good  is  accuracy  if  a  soldier 
doesn't  know  where  to  fire  his  rounds  or  is  so  pinned  down 
by  the  enemy's  superior  fire  power  he  doesn't  have  time  to 
get  a  fix  on  the  target? 

D.  THE  FOOT  SOLDIER  SYSTEM 

The  individual  foot  soldier  is  required  to  both  acquire 
and  engage,  with  small  arms  fire,  enemy  ground  targets. 
Under  target  acquisition  he  is  responsible  for  target 
detection  (the  determination  of  the  existence  or  presence 
of  a  target),  target  identification  (the  determination  of 
the  nature  and  composition  of  the  target),  and  target  loca- 
tion (the  determination  cf  the  coordinates  of  the  target). 
Of  course,  during  the  engagement  phase  he  must  successfully 
neutralize  the  target. 

If  these  combined  processes  could  be  carried  out  with 
zero  error,  every  target  would  be  precisely  located  and 
then  destroyed.   This,  of  course,  is  not  the  case  as  errors 
are  present  in  both  processes  (see  Figure  1,  page  9).   The 
soldier  makes  both  a  target  location  error  and  an  aiming 
error.   There  are  also  errors  in  the  trajectory  of  the  pro- 
jectile due  to  such  variables  as  barrel  wear,  temperature, 
humidity,  and  wind. 

E.  SYSTEM  ERRORS 

This  section  is  an  adaption  to  small  arms  fire  of  errors 
considered  by  C.H.  Hess  [19]  for  artillery  fire.  ■ 
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(XD,YB 


(XA'V 


A    (x^y^) 


«W 


Definition 

target  location  error 

aiming  error 

ballistic  error 

lethality- 


Distribution 


f(XL,YL;x£,y£) 

«*A.Yz;xa,ya) 

f(XB>YB;xb>yb) 


Figure  1.   System  Errors 


1 .  Target  Location  Error 

Targets  to  be  engaged  by  small  arms  fire  are  usually 
located  by  the  foot  soldier  on  the  ground.   These  targets 
are  imprecisely  located  due  to  errors  in  target  acquisition 
and  to  the  fact  that  the  target  is  usually  a  fleeting  target 
The  target's  suspected  location  (XL,YL)  is  distributed  with 
respect  to  the  true  target  location  (x»,y»)  according  to 
some  probability  density  function,  f (XT  ,YT ;x» ,y « ) . 

2.  Aiming  Error 

Because  of  errors  in  firing  techniques  and  correc- 
tions for  other  variables,  an  aiming  error  results  in  a 

separation  of  the  desired  aim  point  (x  ,y  )  and  the  actual 

a   a 

aim  point  (X.,Y.).   The  actual  aim  point  is  distributed 
with  respect  to  the  desired  aim  point  according  to  some 

probability  density  function,  f(X.,Yn;x, ,y  ).   If  there  is 

a   a   a   a 

no  mean  aiming  error,  then  the  expected  aim  point  coincides 
with  the  expected  target  location  (X,  ,Y,  ) . 

3.  Ballistic  Error 

The  impact  of  N  rounds  is  distributed  about  the 
mean  center  of  impact  according  to  some  probability  density 
function,  f  (XR,  YB;x"b  ,y~b  ) . 

4.  Lethality 

The  amount  of  damage  done  to  a  target  by  N-rounds 
is  a  function  of  the  type  target,  the  type  ammunition  fired, 
and  the  distance  from  the  center  of  impact  to  the  target. 
The  lethality  functions  used  in  this  thesis  were  selected 
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based  on  this  knowledge  and  in  all  cases  were  treated  as 
known  analytic  functions. 
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II.   ONE  DIMENSION 

A.  INTRODUCTION 

The  problem  addressed  is  that  of  maximizing  the  proba- 
bility of  hitting  and  thus  killing  a  sharply  defined  one 
dimensional  target.   The  target  was  assumed  to  occupy  a 
sharply  defined  portion  of  a  one  dimensional  coordinate 
system.   A  round  impacting  within  the  target  area  produces 
the  full  effects  of  a  hit,  while  a  round  impacting  outside 
the  target  has  little  or  no  effect  on  the  target. 

The  one  dimension  model  is  presented  to  convey  the 
essential  features  of  the  models  used  and  to  provide  a  solid 
foundation  for  the  more  involved  and  complicated  models 
used  later  in  this  thesis. 

B.  KNOWN  TARGET  LOCATION 

The  single  shot  hit  and  kill  probabilities  were  calcu- 
lated for  fire  delivered  against  a  stationary  target  by 
assuming  that  the  distribution  of  the  impact  points  was 
known  or  could  be  obtained  from  experimental  data.   The 
single  shot  hit  probability,  POQU)  is  obtained  by  use  of 

OOIi 

Equation  (1)  while  Equation  (2)  is  used  to  obtain  the  single 
shot  kill  probability,  P^, 


PSSH     ;   fX  (x)dx  (1) 

~Li  J- 


PSSK  =  ;^(xi)fx  (x)dx  (2) 
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For  both  calculations,  fv  (x)  is  the  density  of  the  impact 
points  while  l(Xj)    is  an  appropriate  lethality  function. 
The  lethality  function  was  used  to  account  for  the  lethality 
of  the  various  types  of  small  arms  ammunition  available. 
1.   Impact  Point 

The  distribution  of  the  impact  points  from  the 
center  of  the  target  is,  in  many  cases,  normal  and  was 
assumed  such  for  this  thesis.   The  Central  Limit  Theorem 
provides  a  theoretical  basis  for  this  empirical  fact,  since 
many  factors  together  cause  a  soldier  to  miss  the  intended 
aim  point.   The  density  used  was  thus 


c-n_i  ? 


l        \°i 

fy  (x)  =  i exp   I  J 

XI       (2tt)^i 

where  uT  is  the  expected  impact  point  and  cT  is  the  measure 
of  dispersion  about  uT. 

2 .   Lethality  Functions 

The  lethality  function,  defined  as  the  Prob  [kill 
target | round  impacts  at  x]  is  introduced  to  provide  a  means 
of  calculating  P.^,   for  all  types  of  small  arms  ammunition 
available,  both  fragmenting  and  non-fragmenting. 

For  non-fragmenting  ammunition  the  zero-one  lethality 
function,  denoted  by  I    (x),  is  used  to  reflect  the  fact  that 
a  hit  is  required  in  order  to  obtain  a  kill.   Thus,  for  a  one 
dimension  target  of  dimension  2L,  the  lethality  function  is 
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IJx)    = 
c 


1;   if  x  e  [-L,L] 

0;   otherwise  . 


Fragmenting  ammunition,  grenades  and  the  like,  does 
not  necessarily  need  to  impact  on  the  target  to  kill  it. 
Fragmenting  ammunition  can  kill  a  fragment  sensitive  target 
by  throwing  shrapnel  on  the  target.   The  damage  done  by  the 
fragmentation  type  round  depends  on  the  target  location  and 
the  point  of  impact  of  the  round.   Figure  2,  page  15,  from 
BRL  report  15^  displays,  rather  vividly,  why  the  negative 
exponential  function  was  chosen  for  this  thesis.   The 
lethality  function  is  therefore 


-C,(x)  =  exp 


-%{|)2 


where  a  is  a  shape  parameter  derived  from  a  fit  to  experi- 
mental lethality  data  which  has  been  independently  determined 
from  fragmentation  tests  or  other  techniques  and  is  called 
the  lethal  radius.   The  variable  x  is  the  distance  from 
the  center  of  impact  to  the  target  center. 
3.   Non-Fragmenting  Model 

Under  the  above  assumptions  the  single  shot  kill 

probability  is  determined  by  evaluation  of  the  integral 

:-yj  2 


SSK 


/  £  (x)  ■ 

(2tt) 


57 


ax 


"l 


(3) 
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Figure   2.   Round  Lethality  Function 
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Use  of  the  zero-one  lethality  function  reduces  this 
evaluation  to  the  interval  [-L,L]  so  that 


L         ,        "* 


aI 


2 


PqqK  =    /   (1)   i- —  exp   I   x  J  dx,  (4) 

SSK    _L         (27r)^ 


finally 


L  -  yT      -L  -  y 


where 


yT         2, 

4>(uT)  =    /  — i-r  e_t  /2  dt 

1      .oo     (2tt  )^ 


is  the  cumulative  distribution  function  for  the  standardized 
normal  distribution. 

^ •   Fragmenting  Model 

For  this  thesis  all  ammunition,  other  than  ball 
ammunition,  is  considered  fragmenting  ammunition.   Therefore, 


PSSK  =    f   lE{x)fX    (x)  dXj  (6) 

_oo  I 


where 


l~(x)    =  prob  [kill  target  |  round  impacts  at  x] 


The  single  shot  kill  probability,  P„K,  is  found  to  be 
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?SSK  "     (a2  +  lx2)h     SXP         '   '     (?) 

5.   Maximizing  Poov- 

An  examination  of  the  conditions  under  which  PCCLr 

is  maximized  leads  to  an  examination  of  the  partial  deriva- 
3  P  3  P 

tives;  ■*— and  ■*-—- [17].   The  results  of  this  examination 

dy         oo 

are  twofold: 

1.  P„n„  is  always  decreased  when  the  expected  impact 
point  is  moved  away  from  the  target  center. 

2.  P„„„  can  be  increased  when  v    >  L,  i.e.  when  the 
expected  impact  point  is  off  the  target.   It  is 
not  unreasonable  to  postulate  that  there  is  no 
net  effect  on  the  expected  impact  point  from 
ballistic  sources.   Therefore  the  condition  that 
u  >  L  must  be  due  to  bias  in  either  target  loca- 
tion, aim  point  selection  or  bot  .       =  ,  1.  . 
can  be  improved  by  increasing  the  dispersion  of 
rounds  about  the  expected  impact  point  when  these 
biases  are  present. 

C.   UNKNOWN  TARGET  LOCATION 

The  results  of  the  last  section  indicate  that  to  effec- 
tively evaluate  weapon  systems,  a  model  which  includes  target 
location  error  is  needed.   Comparing  weapon  systems  based 
on  their  performance  against  stationary,  clearly  visible 
targets  is  somewhat  less  than  satisfactory,  for,  except  in 
rare  instances,  targets  do  not  remain  motionless  and  in 
clear  view. 
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1 .  Target  Location  Error 

The  results  of  the  last  section  merely  reflect  the 
fact  that  in  war,  targets  are  seldom  located  where  they  are 
thought  to  be  located.   The  suspected  location  is  distributed 
with  respect  to  the  true  target  center  according  to  some 
probability  distribution.   For  the  one  dimension  case,  the 
target  location  error  is  assumed  to  be  distributed  normal 
with  mean  u«  and  variance  a»  . 

2 .  Ballistic  Error 

The  impact  of  rounds  about  the  target's  suspected 
location  is  also  assumed  to  be  normally  distributed  so  that 


fX<V 


kB  (2iryab 

where  the  ballistic  dispersion  a,  is  assumed  to  be  composed 
of  terminal  ballistics  dispersion  only.  The  dispersion  due 
to  recoil  was  assumed  to  be  zero. 

3.   Lethality  Function 

The  zero-one  lethality  function  is  used  for  the  non- 
fragmenting  case  while  the  negative  exponential  is  assumed 
for  the  fragmenting  case. 

k.      Analytical  Model 

The  single  shot  hit  and  kill  probabilities  discussed 
earlier,  although  important,  are  not  the  entire  story.   The 
probability  of  hitting  and/or  killing  a  target  when  more  than 
one  round  is  fired  is  now  examined  by  use  of  a  salvo-fire  model 
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In  a  salvo-fire  model  it  is  assumed  that  (1)  the  aim 
point  is  constant  for  all  N-rounds  and  (2)  the  N-rounds  are 
fired  simultaneously  [17].   This  model  describes  numerous 
tactical  situations  among  which  are  (1)  a  soldier  sequen- 
tially firing  N  identical  rounds  at  the  same  aim  point 
(suspected  target  location)  and  (2)  a  squad  of  t  -  men 
firing  N/t  identical  rounds  simultaneously  at  the  same  aim 
point  from  approximately  the  same  location, 
a.   Non-fragmenting  Case 

The  basic  model  and  assumptions  for  salvo- 
firing  of  N-rounds  of  non-fragmenting  ammunition  when  all 
rounds  are  sequentially  fired  at  the  same  aim  point  are: 
(1)  a  soldier  engages  a  target  at  location  x » ,  (2)  he  aims 
his  rifle  at  Xy  where  x«rXy  due  to  target  location  error, 
(3)  he  fires  N-rounds  at  the  suspected  target  location  with 
the  i    round  impacting  at  x,  . ,  and  (4)  the  impact  of  the 
rounds  about  the  fixed  aim  point  are  statistically  indepen- 
dent.  Using  a  rectangular  coordinate  system  (see  Fig.  3) 
with  the  target  reference  point,  x, ,  located  at  0,  the 
assumed  aim  point  is  XT  the  expected  target  location. 

Li 


Lbi  "b2   "L   b3 
Figure  3.   ID  Salvo  Model 
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The  probability  of  hitting  the  target  at  least 
once  with  a  salvo  of  N-rounds  is  calculated  in  the  following 
manner: 


STEP  1.   Compute  the  single  shot  conditional 
probability  of  hitting  the  target 
given  it  is  thought  to  be  at  XT  by 


1        lb 

pssh/xt  =  '  rzsrr  exp    l        J  dxb        (8) 


SSH/XL     (2ir)V 

b 


STEP  2.   Compute  the  conditional  probability 
of  hitting  the  target  at  least  once 
with  a  salvo  of  N-rounds  given  it  is 
thought  to  be  at  XT  by 


(N)  N 

VXL      -  1  -  (1  -  PsSH/XjT  '  (9) 


STEP  3.  Compute  the  unconditional  probabilj 
of  hitting  the  target  at  least  once 
with  a  salvo  of  N-rounds  by 


CO 

PR  (N)  =  1  -    /  (1  -  PSSH/X  )Nfx  (xL)dV  (10) 
-00  JkL     L 


The  probability  of  killing  the  target  with  a 
salvo  of  N-rounds  is  equal  to  P„(N),  since  the  zero-one 
lethality  function  is  used.   The  final  equation  becomes 
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-h 


L   1         I  °h  I     N 
PK(N)=PH(N)  =  1-/  (1-  /  — i—  exp   I   D  J  cb^) 


.  2 


-L   2tt  a, 

D 


(2tt) 7a£ 


dx£  (11) 


for  which  no  closed  form  solution  has  been  found. 
b;   Fragmenting  Case 

The  basic  model  and  assumptions  for  salvo-firing 
of  N -rounds  of  fragmenting  ammunition  when  all  rounds  are 
sequentially  fired  at  the  same  aim  point  are  identical  to 
the  four  of  the  non-fragmenting  case  with  the  additional 
assumption  that  (5)  cumulative  damage  is  negligible. 

The  probability  of  killing  the  target  with  a 
salvo  of  N-rounds  is  determined  as  follows : 


STEP  1.   Con      the  single  shot  conditional 
probability  of  killing  the  target 
given  the  aim  point  is  xL  by, 


oo 


'qcv/v  ~   /    t    exp 


SSK/XT. 


L    -°°  (2tt)20. 
b 


thus 


PSSK/XL  = 
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STEP  2.   Compute  the  conditional  probability 
of  killing  the  target  with  a  salvo 
of  N-rounds  given  the  target  is 
though  to  be  at  XT  by. 


PK/XL(N)  =  1  -  (1  -  P 


ssk/ 


X 


)N . 


(1*0 


L 


STEP  3.   Compute  the  unconditional  probability 
of  killing  the  target  with  a  salvo 
of  N-rounds  by, 


PK(N)  -   /   f 


x£(xL)  PK/xL(N)  dx   , 


(15) 


thus 


N   /  rt         °° 


Under  the  above  assumptions  this  equation  reduces  to 


N 


PK(H) 


/ 


a  X   I,: 
K=l 


1  K-l 


n 


K  +  %2>' 


-Ky 


exp 


2(a2  +  cb2  +  Kc^2) 


(a2  +  ab2  +  Kot2)h 


(17) 
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When  it  is  assumed  that  the  location  error  is 
symmetric  about  the  true  target  location,  i.e.  E(X  )  =  x»  = 

u   =0,  then  the  above  closed  form  equation  reduces  to 

a    ■  * 


N 
P„(N)  =  I 

K       K=l 


-a 


(.*  ♦  o^\ 


K-l 


p       p         p  V, 

(a     +  c^   +  Kc^  ) 2 


(18) 
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III.   TWO  DIMENSION 

A.  INTRODUCTION 

The  addition  of  a  second  dimension  brought  the  model 
much  closer  to  reality  and  provided  a  method  of  comparison 
with  tabulated  data  for  square  targets  [7].   Both  range  and 
deflection  error  were  taken  into  account  by  use  of  bivariate 
error  distributions. 

B.  TARGET  LOCATION  ERROR 

The  target  location  error  was  assumed  to  be  uncorrelated 

bivariate  circular  normal  so  that, 

(XL-x£)2  +  (YL-y^2 
-h   2 

fX   Y  ^l^ n)    =  ~2    exp  °l 

L5  L  ^7Tai> 

C.  BALLISTIC  ERROR 

The  ballistic  error  was  also  assumed  to  be  uncorrelated 
bivariate  circular  normal  so  that, 

(xrV2  +  (^.-V2 


-h  2 

a. 


«     /     \      1  b 

fY  Y  (xh>yh)  =  p  eXP 

b 


D.   AIMING  ERROR 

The  aiming  error  was  assumed  to  be  uncorrelated  bivariate 

2    2    2 

circular  normal  with  mean  (u  ,  u  )  and  variance  o     =o     =o 

x   y  a   x   y 
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E.   LETHALITY 

1.   Non-Fragmenting 


1  ;   if  round  impacts  on  the  target 


0  ;   otherwise 


where  the  target  is  assumed  to  be  of  area  A. 
2 .   Fragmenting 

.  (yxb)2  +  (yL-yb)2 

? 

Vvyb)  =  exp 

F.   ANALYTICAL  MODEL 

1.   Non-Fragmenting  Case 

Using  the  "diffuse"  target  approximation  proposed  by 

J.  Von  Newmann  [4],  and  as  discussed  in  Appendix  A,  the 

conditional  probability  of  killing  the  target  with  one  round, 

given  the  target  location  is  (xT ,yr  )  and  aim  at  (u  ,y  )  is 

j <   J i  X   y 

found  to  be 

,   (mx-*l)2  +  (yyL>2 

2         "^       2      2 
b 

where  the  parameter  W  must  be  adjusted  as  discussed  in 

Appendix  A. 

The  probability  of  killing  the  target  with  N-rounds 

when  the  i    round  is  aimed  at  (y   ,  y   )  is  obtained  from 

xi   yi 
solving  the  integral 
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k.I      yi 


N 


"IT  (1" 
1 


"VV.'  -i   .;   ^XL>YL^        M    ^        iSSK/XL,YL^xi^yi 


xi 
dx£,dy£ 


where  the  origin  of  the  (x,y)  coordinate  system  is  the  center 
of  the  true  target's  location. 

For  salvo -fire,  i.e.  all  N-rounds  are  aimed  at  the 
same  aim  point,  (p^yy)  ,  it  follows  that 


PK(N)    =    PK(yx,yy) 


1    ~    ;      ;    fXr,YT(1    "    PSSK    XT,YT^X>V} 


N 


now 


dx^dy^ 


^    ~   PSSK   XT  ,YT  ^xjVV  ^ 

J-J  Li 


N 


-i* 


1    - 


vr 

ZT~ 2  exP 

W     +  a. 


(yx-xL)2  +   (yy-yL)2 

TT2    .         2 
W      +   a, 
b 


N 


N 

E 

K=0 


-W 


tt2  2 

W     +   a 


K 
K  "2 

exp 


(yx-xL)2  +   (uy-yL)2 

2  2 

w   +  % 
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thus 


PK(N) 


N 

W   E 

K=l 


-W 


2    2 
b 


K-l 


K  (^x-XL)2  +  (VyL)2 
"2 


exp 


W2  +  a.  2  +  K(a  2+a  2) 


W2  +  ob2  +  K(0,2+aa2) 


The  assumption  that  the  constant  aim  point,  (yv3y  ), 

x  y 

selected  by  the  soldier  is  also  the  expected  target  location, 
(Xj  ,Y, ) ,(  reduces  the  above  formula  to 

K-l 


PK(N) 


N 


K=l 


-W 


2      2 
W  +  a, 
b 


W 


W2  +  cb2  +  K(aa2  +  c£2) 


(23) 

Substituting  the  adjusted  value  of  the  parameter  W 
(see  Appendix  A)  into  the  above  formula  reduces  the  formula 
to 


PK(N) 


=   E 
K=l 


-4LV 


klf/2i   +  c. 


K-l 


Hi}/; 


4L2/2tt  +  a  2  +  K(o  ^+C02) 
b       a   £ 


(24) 


for  a  square  target  of  sides  2L.   The  computer  program  used 
to  calculate  PK(N)  for  N=5,10  (both  with  and  without  the 
adjustment  on  the  parameter  W)  is  given  in  Appendix  B. 

An  indication  of  the  nature  of  the  diffuse  target  approxi- 
mation was  obtained  by  comparing  the  computer  calculated  Pv(N) 
values  with  the  exact  values  for  a  square  target  of  sides  2L. 
The  latter  values  are  tabulated  in  Reference  7.   Figures  k 
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and  5  graphically  depict  this  comparison  for  N  =  5,  10,  and 

W  adjusted  as  per  Appendix  A.   Table  1  is  given  to  reflect 

the  improvement  of  the  approximation  after  adjustment  of  the 

parameter  W.         0   0   ,  . 

[2(aa2+a£2)]^  ab  (2)* 

The  ratios   and  s used  in  Figures 

L  L 

4  and  5  were  formed  in  order  that  a  direct  comparison  could 

be  made  with  the  tabulated  results  of  Reference  7.   These 

ratios  were  formed  by  multiplying  through  equation  (24)  by 

((2)Vl)2  to  obtain 


N  M 

PK(N)  =   £   £ 
K       K=1  W 


-4/tt 


4A  +  [(^!b32 


L 


L 


K-l 


4/fT 


Li  Li 


(25) 


The  presence  of  the  (2)2  factor  in  the  multiplier, 


U   2 

(2)fL)  ,  was  found  to  be  necessary  in  order  to  obtain  proba- 
bilities close  to  those  tabulated  in  Reference  7  and  noted 
by  H.  J.  Helgert  in  Reference  6.   Although  there  was  no 
a  priori  reason  to  include  the  (2)  factor  in  these  ratios, 
calculations  which  did  not  include  it  in  the  ratios  resulted 
in  salvo  kill  probabilities  far  below  those  of  Reference  7. 
2,   Fragmenting  Case 

The  use  of  the  negative  exponential  lethality 
function  for  fragmentation  ammunition  resulted  in  the  final 
formula 
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0.9 


FIGURE  h.    EXACT  vs  DIFFUSE   (N=5) 

Exact  (7) 

Diffuse 
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0.9 


P(N) 

^.8J 


FIGURE  5.  EXACT  vs  DIFFUSE   (N=10) 


Exact  (7) 


0.7  - 


0.6 


0.5 


0.4 


0.3_ 


0.2 


0.1 


Diffuse 


L 
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NO  ADJUSTMENT 
RATA     RATB 


VL 


c.W)* 


PK(N) 


RATA 


(2)'2a 


ADJUSTED 

RATB     PK^N) 
C2(aa2+a/)]^ 


TABULATED 
PK(N) 


11 


16 


1 

0.8715 

2 

0.6902 

4 

0.3072 

7 

0.1187 

11 

0.0506 

16 

0.0244 

1 

0.5055 

2 

0.4974 

4 

0.2674 

7 

0.1124 

11 

0.0494 

16 

0.0241 

1 

0.1774 

2 

0.2239 

li 

0.17':0 

7 

0.0925 

11 

0.0452 

16 

0.0230 

1 

0.0633 

2' 

0.0882 

4 

0.0860 

7 

0.0620 

11 

0.0365 

16 

0.0206 

1 

0.0263 

2 

0.0379 

4 

0.0422 

7 

0.0360 

11 

0.0258 

16 

0.0167 

1 

0.0125 

2 

0.0183 

4 

0.0213 

7 

0.0201 

11 

0.0166 

16 

0.0123 

11 


16 


1 

O.9631 

2 

.8872 

il 

.5138 

7 

12222 

11 

.0985 

16 

.0482 

1 

.6796 

2 

•  7059 

4 

.4567 

7 

.2111 

11 

.0963 

16 

.0476 

1 

.2732 

2 

•  3510 

4 

•  -J      J 

7 

.1754 

11 

.0882 

16 

.0456 

1 

.1017 

2 

.1441 

4 

.  600 

7 

.1189 

11 

.0716 

16 

.0407 

1 

.0429 

2 

.0628 

4 

.0776 

7 

.0695 

11 

.0507 

16 

.0330 

1 

.0205 

2 

.0305 

4 

.0395 

7 

.0391 

11 

.0327 

16 

.0244 

0.9611 

.9804 

.5257 

.2245 

.0989 

.0481 

.6337 

•  7139 

.4659 

.2131 

.0966 

.0476 

.2359 

.3443 

.  3120 

.1767 

.0885 

.0456 

.0854 

.1389 

.H 

.1195 

.0718 

.0407 

.0357 

.0602 

.0776 

.0697 

.0508 

.0330 

.0170 

.0292 

.0394 

.0391 

.0327 

.0244 

TABLE  I 
SALVO  KILL  PROBABILITIES  FOR  N=10 
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PK(N) 


2N 
K=l 


K 


-a 


2  j.  ~   2 
a  +  a, 

D 


K-l 


K  (VXL)2  +  (VyL)2 
2  .2  .  _  2  ,  „,     2,   2 


exp 


■  +  ab  +  K(aa  +a£  ) 


a2  +  ab2  +  K(aa2+a/) 


,  (26) 


where  a  is  the  lethal  radius  of  the  round.   The  similarity 
between  this  final  equation  and  that  of  the  non-fragmenting 
case  were  expected.   The  diffuse  target  approximation  used 
a  negative  exponential  function  similar  to  the  lethality 
function  used  in  this  case.  The  parameter  a  in  this  case  is  a 
function  of  both  the  type  ammunition  used  and  the  type  target 
engaged,  while  the  parameter  W  in  the  diffuse  target 
approximation  is  a  function  of  the  target  and  not  of  the  type 
ammunition  fired. 
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IV.   OPTIMAL  AIMING 

When  firing  at  a  target  whose  true  location  is  unkown , 
but  a  "guess"  can  be  made  as  to  the  suspected  target  loca- 
tion the  question  of  where  to  aim  arises.  For  the  case  when 
two  rounds  are  to  be  fired  and  assuming  that  (1)  the  origin 
of  the  coordinate  system  is  at  the  center  of  the  suspected 
target  location  and  (2)  that  the  range  error  is  of  greater 
importance  than  deflection  error,  the  optimal  two  aim  points 


are  u   ,  and  u   where  u 


y 


v 


,  as  shown  by  Figure  6. 


Figure  6 


The  probability  of  killing  the  target  is 


-k 


PK(y,-u) 


2W 


„2,    2      2      2 
W  +  on      +  a       +   a. 
I  a     b 


-y 


..2      2~         2      2 
W  +  c    +  an      +   a. 
exp  a     I  b 


(27) 


-V 


^  2    2 

_wf W  +a, 

2    2    2    2   ®XP  , 


(W2+cb2)[2(oa"+a^)+W"+cb" 


v;here  y  .   =  -y   =  y  . 
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For  the  case  when  three  rounds  are  to  be  fired  the  optimal 


aim  points  are  \i      ,  0,  u   as  depicted  in  Figure  7 

yl      y2   v 


y-!  /  y- 


^2 
Figure  7 


The  probability  of  killing  the  target  is  then 


PR(y  ,0,-v0 


W2 


2\ 


£ 


lA-a  2+a  2+a  2      w+a^+a  2+a  2 
a      I      b  £      a      b 


,  f  -„2    1 

ir+a  2+a01 
exp      t       a      J?,j 


(  ,,2 


-\ 


o        o        p  p 

j(cf>.       -2)      u2 

-"2 


VJ  exp 


W2-^  2 


(W2+cb2)[2(a/+oa2)+W2+0b2] 


/2+cb2)(2c  2+2o/+obW: 
2W  exp 

(W2^2 )  [2 (o  2-kJt2  HW2-^2 ] 


w2 


W  exp 


(28) 


(W^2 )  [2  (aa2+c£2  )+VT+ab2 


3^ 


Under  the  above  assumptions,  the  case  when  IT  rounds  are 
to  be  fired  at  u    and  K  rounds  at  u   ,  where  N  >  2 .  was 

y1  y2       = 

not  pursued  due  to  the  problems  of  integrating  the  product 
of  two  sums  over  a  double  integral.   The  equation  to  be 
solved  is 


00  00 

a  f3L,Y7 

-00  -CO     Tj5   L 


N 

E 

K=0 


-w2 


i  t2-l  2A      2 
£   a 


K 


N 
X  E 

K=0 


-W2 


W  +0.  +a 
A   a 


I 


k  (VJY-yL)2+xL2 

~2   2   2   2 
exp       £   a 


K 


^Y-yL)+xL 


«Lal 


K     2  WW2 
exp       £   a 


dx^ 


(29) 


Thus,  when  the  true  target  location  is  unknown,  in 
attempting  to  maximize  the  probability  of  hitting  the  tarL 
and  thus  of  incapacitating  it,  the  rounds  should  be  delivered 
so  as  to  cover  both  the  suspected  target  center  and  the  area 
immediately  surrounding  the  target.   This  can  be  accomplished 
at  short  ranges  by  using  grenade,   shotgun  or  flechette  type 
ammunition,  while  at  longer  ranges  ball  ammunition  fired  from 
a  light  machine  gun  or  similar  rapid  fire  system  would  be 
appropriate. 
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V.   CONCLUSIONS 

When  the  target's  true  location  is  known ,  i.e.,  no 
target  location  error,  the  optimal  aim  point  is  one  which 
places  the  expected  impact  point  of  the  round  on  the  target 
center  (see  Maximizing  ?„„„).   Therefore,  as  is  obvious, 
a  desirable  weapon  system  is  one  which  places  a  round  on 
the  target  center  with  a  probability  near  one.   A  high  degree 
of  accuracy  and  a  slow  rate  of  fire  are  two  characteristics 
of  such  a  weapon  system. 

However,  except  in  very  rare  instances,  targets  do  not 
remain  motionless  and  in  clear  view.   More  often  than  not, 
the  target  appears  as  a  fleeting  target  whose  true  location 
is  unknown.   Under  these  more  realistic  conditions  dispersion 
about  the  suspected  target  location  is  desirable.   Depending 
on  the  range  from  the  firer  to  the  target  this  can  be 
accomplished  by: 

(a)  Spraying  the  target  area  with  ball  ammunition  from 
a  fairly  inaccurate  but  rapid  firing  system. 

(b)  Firing  fragmentation  ammunition,  whether  it  be 
shotgun,  grenade,  or  flechette,  at  the  suspected  target 
location. 
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APPENDIX  A 
DIFFUSE  TARGET  APPROXIMATION 

To  obtain  a  solution  in  closed-form  for  a  sharply  defined 
target  the  diffuse  target  approximation  is  employed,  i.e. 
rounds  impacting  near  the  target  produce  a  nearly  full  effect 
on  the  target,  while  those  falling  far  from  the  target 
produce  a  very  small  effect  on  the  target.   This  approach  was 
originally  suggested  by  J.  von  Neumann  [4]. 

The  negative  exponential  function  Is  used  for  this 
approximation.   Thus, 

(x,-x,  )2  +  (yr-yj2 


JL^b7   '  ^L  Jb 

W2 


fX,,Y  ^VV  =  eXP 
d'  d 


The  negative  exponential  is  used  since  (1)  th     ;ative 
exponential  drops  off  sharply  outside  the  target  area  and 
(2)  it  greatly  simplifies  the  calculations. 
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Figure  A-l  shows  the  original  situation  while  Figure  A-2 
shows  the  approximation. 


Figure  A-l 


2L 

Y 

V 

■L 

xc. 

+T. 

Figure  A-2 

The  diffuse  target  approximation  of  Figure  A-2  allows 
the  integration  to  be  carried  out  from  -00  to  +00  with  a 
minimum  amount  of  probability  in  the  tails,  i.e.  beyond 
Xp+L  and  X.-L.   To  further  minimize  the  discrepancies 
between  the  original  situation  and  the  approximation  the 
parameter  W  must  be  adjusted. 
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Adjusting  W  so  that  the  areas  of  Figure  A-l  and  A-2 
are  equal  implies 

(xrxb)2  +  (y£-yb)2 

r       r  w  p 

J    J   exp  dxdy  =  4L    . 

—  00    —  oo 


Performing  the  above  integration  implies  that 


w2 

2lt 

= 

kL2 

w2 

= 

2 

27T 

w2 

Target    Area 

277 
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APPENDIX  B 
COMPUTER  OUTPUT  FOR  A  SALVO  OP  5  ROUNDS 

This  table  compares  the  results  of  the  diffuse  target 
approximation  for  a  salvo  of  N  =  5  rounds. 
The  inputs  are; 

N  =  NRDS  =  Number  of  rounds  in  the  salvo 


RATA      =   gj  °b 


[2<°a2  +  °£2if 


RATB 


BINOMC(N)=   Binomial  coefficients 


See  Appendix  D  for  the  computer  program. 
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PK(N) 

N 

RATB 

RATA 

DIFFUSE 

EXACT 

5 

1. 

1. 

0.872415 

0.893900 

5 

1. 

2. 

0.67255^ 

0.709200 

5 

1. 

4. 

0.302942 

0.3H400 

5 

1. 

7. 

0.118124 

0.119300 

5 

1. 

11. 

0.050578 

0.050700 

5 

1. 

16. 

0.024403 

0.024300 

5 

2. 

1. 

0.544349 

0.525200 

5 

2. 

2. 

0.493774 

0.510400 

5 

2. 

4. 

0.264370 

0.270800 

5 

2. 

7. 

0.111875 

0.113000 

5 

2. 

11. 

0.049400 

0.049500 

5 

2. 

16. 

0.024125 

0.024100 

5 

4. 

1. 

0.204383 

0.184500 

5 

4. 

2. 

0.229438 

0.229500 

5 

4. 

4. 

0.173284 

0.175800 

5 

i». 

7. 

0.092206 

0.092900 

5 

4. 

11. 

0.045181 

0.045300 

5 

4. 

16. 

0.C23075 

0.023000 

5 

7. 

1. 

0.074655 

O.O658OO 

5 

7. 

2. 

0.091964 

0.090400 

5 

7. 

4. 

0.088312 

0.038800 

5 

7. 

7. 

O.O61983 

0.062300 

5 

7. 

11. 

0.036566 

0.036600 

5 

7. 

16. 

0.020606 

0.020600 

5 

11. 

1. 

0.031285 

0.027300 

5 

■11. 

2. 

0.039817 

0.038900 

5 

11. 

ii , 

0  .  0  4  2  r:  9  0 

0.042600 

5 

11. 

7- 

0.02: 

0 . 0  3  c  1  c : 

5 

11. 

11 . 

0.025803 

0.025800 

5 

11. 

16. 

0.016701 

0.016700 

5 

16. 

1. 

0.014973 

0.013000 

5 

16. 

2. 

0.019296 

0.018800 

5 

16. 

4. 

0.021605 

0.021500 

5 

16. 

7. 

0.020219 

0.020200 

5 

16. 

11. 

0.016619 

0.016600 

5 

16. 

16. 

0.012319 

0.012300 
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APPENDIX  C 
COMPUTER  OUTPUT  FOR  A  SALVO  OF  10  ROUNDS 

This  table  compares  the  results  of  the  diffuse  target 
approximation  for  a  salvo  of  10  rounds  with  the  tabulated 
values  of  Reference  7. 

The  inputs  are: 

N  =  NRDS  =  Number  of  rounds  in  the  salvo 

h 

(2)K 


RATA     =    L 
RATB      = 


[2  (a2  +  ?  2] 


BIN0MC(N)=  Binomial  coefficient! 


See  Appendix  D  for  the  computer  program 
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PK(N) 

N 

RATB 

RATA 

DIFFUSE 

EXACT 

10 

1. 

1. 

0.963104 

0.961100 

10 

1. 

2. 

0.887269 

0.908400 

10 

1. 

4. 

0.513899 

0.525700 

10 

1. 

7. 

0.222290 

0.224500 

10 

1. 

11. 

0.098597 

0.098900 

10 

1. 

16. 

0.048210 

0.048100 

10 

2. 

1. 

0.679618 

0.633700 

10 

2. 

2. 

0.705957 

0.713900 

10 

2. 

4. 

0.456730 

0.465900 

10 

2. 

7. 

0.211172 

0.213100 

10 

2. 

11. 

0.096356 

0.096600 

10 

2. 

16. 

0.047668 

0.047600 

10 

4. 

1. 

0.273233 

0.235900 

10 

4. 

2. 

0.351C53 

0.344300 

10 

4. 

4. 

0.308397 

0.312000 

10 

4. 

7. 

0.175413 

0.176700 

10 

4. 

11. 

0.088294 

0.088500 

10 

4. 

16. 

0 . 045616 

0.045600 

10 

7. 

1. 

0.101787 

0.085400 

10 

7. 

2. 

0.144124 

0.138900 

10 

7. 

4. 

0.160043 

c. 160600 

10 

7. 

7. 

0.118929 

0.119500 

10 

7. 

11. 

0.071679 

0.071800 

10 

7. 

16. 

0.040776 

0.040700 

10 

11. 

1. 

0.042912 

0.035700 

10 

11. 

2. 

0.062891 

0.060200 

10 

11. 

h 

0.07'  7    " 

0.077600 

10 

11 . 

7. 

0.069592 

0.069  700 

10 

11. 

11. 

0.050724 

0.050800 

10 

11. 

16. 

0.033091 

0.033000 

10 

16. 

1. 

0.020583 

0.017000 

10 

16. 

2. 

0.030568 

0.029200 

10 

16. 

4. 

0.029568 

0.039400 

10 

16. 

7. 

0.039100 

0.039100 

10 

16. 

11. 

0.032731 

0.032700 

10 

16. 

16. 

0.024437 

0.024400 
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COMPUTER  PROGRAM  FOR  TWO  DIMENSION  MODEL 
WITH  CIRCULAR  NORMAL  DISTRIBUTIONS 


C 

c 
c 

c 
c 


1300 

1002 
1003 


100 


101 


104 

1001 

103 

99 


THI 
THE 

DIS 

IMPLICI 

THE 

NUM 

DIMENSI 

READ(5, 

READ( 5« 

READ(5, 

READI5, 

FORMAT 

FORMAT 

FORMAT 

NRDS=5 

N=NRDS 

DO  99  I 

DO  103 

DO  100 

TFRM2A= 

TERM2=( 

TERM3A= 

TFRM3=1 

TERMCK) 

CONTINU 

SUM=0.0 

PL  101 

SUM=SUM 

CONTINU 

P  K I  L  L  =  S 

GO  TO  1 

W  R  1 1  E  I  0 

FORMAT ( 

CONTINU 

CCNT  INU 

STOP 


S  PROGRAM  CCMPUTES  SALVC  KILL  PROBABILITIES  FOR 
TWO  DIMENSION  SALVO  MODEL  WITH  CIRCULAR  NORMAL 

TRIBUTI ONS 

T  REAL«8(  A-H,C-Z) 
FOLLOWING  MUST  BE  DIMENSIONED  TO  NRDS.THE 

BER  CF  RCUNDS  FIRED  IN  THE  SALVO 

ON  TER^K  5) , RATA (6)  , RATE (6, 6)  , BI NOMC ( 5 ) ,T AB ( 6 , 6 ) 

1000  )  ( RATA(J)  ,  J  =  l, 6}  ' 

1JJD)  (  (RAT3(y,L)  ,L=1,6) 

1002)  (  BINOMC(K)  ,K=1,5) 

1003)  ( (TAB(M,L) ,L=1,6) t 
( 18F4.0) 
( 15F6.0) 
(8F9.6) 


,M=1,6) 
1=lf  6) 


=  1. 
J  =  l 

K  =  l 
1.2 
-1  . 
TER 
.27 
=  BI 


,6 

,N 

73  +  ( RATB(  I, J  J  )**2 

273/TERM2A)**(K-l) 

M2A+(K'-MRATA(  I) )**2J 

3/TERM3A 

NOMC (K)*TERM2*TERM3 


K = 1  «  *  J 

+  TERM<  K) 
c 

UM 

34 

6,1001)  NRDS,RATA(I  )  - 

2  IX 

E 


{  1  ,  J)  ,PKILL,TA 

I  3 , 5X , F4 . 0 , 5X «  F4 . 0 , 5X » F9 . 6 , 5X , F  9 . 6 ) 


(I  ,  J) 
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